Evidence for a new light spin-zero boson from cosmological gamma-ray propagation? 
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Recent findings by imaging atmospheric Cherenkov telescopes indicate a large transparency of the 
Universe to gamma rays, which can be hardly explained within the current models of extragalactic 
background light. We show that the observed transparency is naturally produced by an oscillation 
mechanism — which can occur inside intergalactic magnetic fields — whereby a photon can become a 
new spin-zero boson with mass m >C 10 -10 eV. Because the latter particle travels unimpeded 
throughout the Universe, photons can reach the observer even if the distance from the source 
considerably exceeds their mean free path. We compute the expected flux of gamma rays from 
blazar 3C279 at different energies. Our predictions can be tested in the near future by the 
gamma-ray telescopes H.E.S.S., MAGIC, CANGAROO and VERITAS. Moreover, our result 
provides an important observational test for models of dark energy wherein quintessence is coupled 
to the photon through an effective dimension-five operator. 
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Thanks to the advent of imaging atmospheric 
Cherenkov telescopes (IACTs) like H.E.S.S., MAGIC, 
CANGAROO III and VERITAS, the detection of pho- 
tons from astrophysical sources in the Very High En- 
ergy (VHE) window above a threshold close to 100 GeV 
and up to some lOTeV has become an exciting real- 
ity. Active galactic nuclei (AGN) are known emitters 
of photons in that energy range through specific pro- 
cesses arising from their accretion [l| . A number of such 
emitters have been observed in the last few years, and 
about 20 AGN at VHE are known today @. The red- 
shift z at which thesesources are detected is steadily in- 
creasing, and recently the MAGIC collaboration has re- 
ported [3j an impressive record at z = 0.538 with the 
discovery of the emission from 3C279 0]. 

According to the current understanding, the flux of 
VHE gamma rays from distant sources is attenuated 
in an energy-dependent way by the interaction with 
background photons in the Universe [f|. Specifically, 
the dominant process for the absorption of a photon 
of energy E is pair-production 77 — > e + e~. Actually, 
(7(77 — > e + e _ ) becomes maximal for a background pho- 
ton energy e ~ (500 GeV/E) eV. In the energy interval 
100 GeV < E < 1 TeV the absorption is dominated by 
the interaction with optical/infrared photons of the meta- 
galactic radiation field (MRF), often called extragalactic 
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background light (EBL) 6], while the interaction with 
the CMB becomes dominant at E ~ 10 6 GeV. 

A general trend seems to emerge from the most recent 
determinations of the EBL density @, H[ • Observations 
of AGN at redshift z = 0.165 and z = 0.186 entail that 
intergalactic space might be more transparent to gamma 
rays than previously thought. Since only basic QED and 
the theory of special relativity are involved in the calcula- 
tion of 17(77 — > e + e~), it might be difficult to account for 
the observed attenuation without invoking new physics 
beyond the Standard Model [9j]. Even in a model with 
the lowest EBL density compatible with cosmology, one 
expects for a source at redshift z ~ 0.538 an attenuation 
by less than 0.5 at 100 GeV, which increases to about 
e~ 4 ~ 0.018 at E = 500 GeV Q5- This makes a source 
like 3C279 hardly visible in the VHE range — or, in any 
case, causes a severe difference in the attenuation of the 
signal at, say, 100 GeV, as compared to the attenuation 
at, say, 500 GeV. Yet, MAGIC has found that the sig- 
nal from 3C279 collected in the region E < 220 GeV has 
more or less the same statistical significance as the signal 
in the range 220 GeV < E < 600 GeV (6.1<r in the former 
case, 5.1cr in the latter) Q. It seems therefore interesting 
to explore which kind of new physics would be required 
to explain the observations of VHE sources. 

In this paper, we show that a way out of the above dif- 
ficulty is provided by a simple mechanism, according to 
which part of the photons are converted into light spin- 
zero bosons, which can travel unimpeded through cos- 
mological distances; part of such light bosons are in turn 
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reconverted into photons and detected, so that every- 
thing goes as if an anomalously small attenuation were 
operative. As a result, photons can reach the observer 
even if their mean free path is considerably smaller than 
the distance from the source. We show that this mech- 
anism indeed predicts an effect which is quantitatively 
adequate to explain the data. 

A natural possibility of this kind arises in extensions 
of the standard model containing a new light spin-zero 
boson X coupled to the photon through an effective 
dimension-five operator. In order to be specific, we shall 
be concerned throughout with the effective Lagrangian 
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where <j) stands for the X field and M is a parameter 
with the dimension of an energy, playing the role of the 
inverse of a coupling constant [ill ]. A Lagrangian like 
C^y appears in a wide class of realistic four-dimensional 
models fl2j , in particular in the phenomenology of axion- 
like particles, as well as in compactified Kaluza-Klein 
theories [l3[ and superstring theories 0|. Moreover, 
it has been argued that the presence of should be 
a generic feature of quintessential models of dark en- 
ergy 15 1 . A particular case concerns the axion 



which is characterized by the mass-coupling relation 
77i - (10 10 GeV/M)eV. 

Before proceeding further, we recall that astrophysics 
provides a lower bound on M, depending on the mass m 
of the X boson. Failure to detect X bosons from the 
Sun by the CAST experiment at CERN sets the robust 
bound M > 1.14 • 10 10 GeV for m < 0.02eV [l7|, which 
practically coincides with the theoretical bound derived 
from the properties of globular clusters [Hj|. Further- 
more, the stronger bound M > 3 ■ 10 11 GeV holds for 
X bosons with m < 10 -10 eV, based both on the ener- 
getics of the supernova 1987a and on observations of 
time-lag between opposite-polarization modes in pulsar 
radio emission [20j . 

Owing to £^ 7 , the interaction eigenstates differ from 
the propagation eigenstates in the presence of a magnetic 
field B, so that j-X interconversion occurs. Coherent 
j-X mixing can be understood as an oscillation process 
quite similar to that taking place for massive neutrinos 
of different flavours, apart from the fact that in the case 
described in this paper the external B field is necessary 
to account for the spin mismatch [2l| . 

As an illustrative example, we compute the fraction of 
the initial photon flux which survives the distance from 
the source 3C279 at z = 0.538 when both the absorption 
from EBL and j-X oscillations in intergalactic magnetic 
fields are taken into account [22| . 

Intergalactic magnetic fields have a complicated and 
poorly known morphology, which reflects the pattern of 
baryonic structure formation as well as its subsequent 
evolutionary history [23| . However, for the present needs 
it is sufficient to suppose that B is constant over a do- 
main of size Ldom, with B randomly changing its direc- 
tion from one domain to another but keeping the same 



strenght. Values to be used throughout are B ~ 10~ 9 G 
and idom — 1 Mpc, which are close to existing upper lim- 
its but consistent with them [24 1. 

Let us consider first the propagation of a photon beam 
over a single magnetic domain, in the presence of a 
cold intergalactic plasma with plasma frequency U3 V \ = 
10~ n ^/n e /cm~ 3 eV, where n e 
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^4iran e /m e ~ 3.69 
notes the electron density. In order to have an unsup 
pressed amplitude, we work in the strong-mixing regime 
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pi 



\M/2B. A conservative 



estimate of the density of intergalactic plasma yields 
77, e ~ 10~ 7 cm -3 [2||, resulting in the plasma frequency 
w p i ~ 1.17 • 10~ 14 eV [26j]. Therefore, the latter condition 
takes the explicit form Kto/IO" 10 eV) 2 - 1.37 • 10 8 | <C 
0.38(uVGeV)(-B T /10- 9 G) (10 10 GeV/Af). Since we are 
interested in the energy range E > 10 2 GeV, we find that 
for e.g. M > 4 • 10 11 GeV the bound m < 10~ 10 cV has 
to be satisfied [27']. We remark that the present mech- 
anism works for arbitrarily small values of to, provided 
M is considerably smaller than the Planck mass Mp ~ 
1.22 ■ 10 19 GeV 28]. As a consequence, our result also 
applies to models of dark energy wherein quintessence 
enjoys a photon coupling described by [la], thus ul- 
timately providing an important observational test for 
these models. 

It is straightforward to solve the beam propagation 
equation once photon absorption by the EBL is taken 
into account and produces a finite photon mean free path 
A 7 [291 ]. In the strong-mixing regime, the probability for 
a photon to become a X boson after a distance y reads 



(2) 
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FIG. 1: The two lowest lines give the fraction of photons sur- 
viving from a source at the same distance of 3C279 without 
the oscillation mechanism, for the "best-fit model" of EBL 
(dashed line) and for the minimum EBL density compati- 
ble with cosmology [Toj ] . The solid line represents the pre- 
diction of the oscillation mechanism for B ~ 1CP 9 G and 
Z/dom — 1 Mpc within the "best-fit model" of EBL. The gray 
band is the envelope of the results obtained by independently 
changing B and Ldom within a factor of 10 about their pre- 
ferred values. 
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whereas the probability that a photon remains a photon 
reads 
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where we have introduced the dimcnsionlcss parameter 
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The Fazio-Stecker relation between the photon energy 
and the redshift of its horizon [3(| yields A 7 ~ 450 Mpc 
at E — 500 GeV, so that 5 ~ 12.4 for this case study. 

Over distances y 3> Tdom, the transition probabilities 
P-y^a(y) and P 7 ^ 7 (y) arise as the incoherent average of 
Eqs. ([2|) and ((3]) over N ~ (y/Ldom) domains crossed by 
the beam, respectively. Assuming (as before) that the 
beam propagates along the y direction and choosing the 
x and z directions arbitrarily in the orthogonal plane, the 
problem becomes truly three-dimensional, because of the 



random orientation of the magnetic field. Consequently, 
the beam state is described by the vector (p/ x , 7 Z , a) . 

We have written down the propagation equations de- 
scribing the absorption of photons due to the interaction 
with the EBL and their oscillations into X bosons (and 
vice- versa). Similarly to Ref. (29|, we are led to the trans- 
fer equation 
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and is the angle between the x axis and the extra- 
galactic B in a single domain. Starting with an unpolar- 
ized photon beam, we propagate it by iterating Eq. ([5]) 
as many times as the number of domains crossed by 
the beam, taking each time a random value for the an- 
gle 9 (this reflects the random orientation of B), for 
E = 500 GeV and z = 0.538, corresponding to blazar 
3C279. We next repeat such a procedure 10' 000 times. 
Upon averaging over all these realizations of the propa- 
gation process, we find that about 13% of the photons 
arrive to the Earth, resulting in an enhancement by a 
factor of about 20 with respect to the flux expected in 
the absence of the proposed oscillation mechanism; the 
comparison is made with the "best-fit model" described 
in in the first paper in Ref. [13(a)]. The same calculation 



gives a fraction of 76% at 100 GeV (to be compared to 
67% without the oscillation mechanism) and a fraction 
of 3.4% at 1 TeV (to be compared to 0.0045% without 
the oscillation mechanism) (cfr. Fig. |T|) . In addition, we 
have checked the stability of our result against indepen- 
dent variations of B and Tdom within a factor of 10 about 
their preferred values. The resulting spectrum is repre- 
sented by the gray band in Fig. 1. We remark that the 
standard deviation of the above averaging procedure lies 
well inside the gray band. Our predictions can be tested 
in the near future with IACTs. 
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